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Abstract 

We investigate the regulators of elliptic curves with rank 1 in some 
families of quadratic twists of a fixed elliptic curve. In particular, we 
formulate some conjectures on the average size of these regulators. 
We also describe an efficient algorithm to compute explicitly some 
of the invariants of a rank one quadratic twist of an elliptic curve 
(regulator, order of the Tate-Shafarevich group, etc.) and we discuss 
the numerical data that we obtain and compare it with our predictions. 

1 Introduction and notations 

We study the regulators of elliptic curves of rank 1 in a family of quadratic 
twists of a fixed elliptic curve E denned over Q. Methods coming from 
Random Matrix Theory, as developed in |K-S] . |CKRS] . |CFKRS] . etc., allow 
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us to derive precise conjectures for the moments of those regulators. Our 
hope is that these moments will help to make predictions for the number 
of curves with extra-rank (i.e. the number of even quadratic twista^j] with a 
Mordell-Weil rank > 2, or the number of odd quadratic twists with Mordell- 
Weil rank > 3). Then, we describe an efficient method, using Heegner- 
point construction, for computing the regulator (and the order of the Tate- 
Shafarevich group) of an elliptic curve of rank 1 in a family of quadratic 
twists. Finally, we discuss and compare our extensive numerical data (for 
some families of odd quadratic twists of the curves llal, 14al, 15al and 
17al) with our predictions. 

From a numerical and experimental point of view, the situation of odd 
quadratic twists really differs from the one of even quadratic twists. Indeed, 
in the latter case, for each curve E d in a family (E d ) d of even quadratic twists 
of a fixed elliptic curve E, one has to compute the special value L(E d , 1) of 
its L-function at s = 1 and determine if it is zero or not. If L(E d , 1) = then 
the curve E d has extra-rank. Otherwise the curve has rank 0, the regulator is 
simply 1, and the Birch and Swinnerton-Dyer conjecture allows us to deduce 
the value of |lIl(L^)| from that of L(E d ,l). The computation of L(E d ,l) 
is done via a Waldspurger's formula which, roughly speaking, states that 
L(E d , 1) is, up to a fudge factor, the square of the |rf|-th coefficient of a 
weight 3/2 modular form given by an explicit linear combination of theta 
series. It follows that, in this case, computations are possible for very large 
families of quadratic twists (see for example [Rub], |Qua| , etc.). Note that 
the numerical data coming from these computations are in close agreement 
with the well-known conjectures of |CKRS] about extra-vanishing (coming 
from the models of Random Matrix Theory), or on the behavior of the Tate- 
Shafarevich groups Ul(E d ) of E d (see [Qua], [Del] ). 

In the rank 1 case, numerical investigation appears to be much more com- 
plicated and, as far as we know, has never been done before. In that case, we 
first have to compute the value of the derivative L'(E d , 1) for each curve E d in 
the family of odd quadratic twists. However, there is no Waldspurger's for- 
mula to compute this value directly, and furthermore from this value one can 
only deduce (assuming it is non-zero and under the Birch and Swinnerton- 
Dyer conjecture) the value of the product R(E d ) 1111(^)1 where R(E d ) is the 

*An odd (resp. even) quadratic twist of E is a quadratic twist such that the sign of the 
functional equation of its L-function is —1 (resp. +1). By the Birch and Swinnerton-Dyer 
conjecture this is equivalent to say that its Mordcll-Wcil rank is odd (resp. even) 
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regulator of E d . Thus we also need to be able to evaluate at least one of the 
two terms of this product]^] The only (known) efficient way to do this is to 
write down a generator Gd of Ea(Q) and to compute R(Ed) = h{Gd) where 
h is the canonical heighlnof Ed- 



The method we used in this paper is to first adapt the Heegner-point 
construction to our situation in order to construct a generator Gd and then 
replace the Waldspurger's formula by the formula of Gross and Zagier. This 
allows us to compute directly the regulator R(Ed) and at the same time 
the order of the Tate-Shafarevich group |lll(i^)| (assuming the Birch and 
Swinnerton-Dyer conjecture). 

Hypothesis. From now on, we assume the truth of the Birch and 
Swinnerton-Dyer conjecture. 

We now give some notations. Fix an elliptic curve E defined over Q and 
let N be its conductor. The L-function of E is 



L(E, s) = J2 a{n)n~ s , R(s) > 3/2 



n>l 



It is now a classical and deep result that L(E, s) can be analytically continued 
to the whole complex plane and satisfies a functional equation: 



A(E,s):=\^-j T(s)L(E,s)=wA(E,2-s) 

where w = ±1 gives the parity of the order of vanishing of L(E, s) at s — 1. 
Let d be a fundamental discriminant. We denote by Ed the quadratic twist 
of E by d. The curves E and Ed are isomorphic over the quadratic field 
but not over Q. We denote by ipd (V* if d is clear in the context) the 



2 For some families of elliptic curves (Fj)j, there exists a generic point in the Mordell- 
Weil group Fj(Q), thus one can separate the terms in this product and a direct investiga- 
tion is possible (see |De-Du] ) . However, such families for which we know in advance the 
regulator are very special and in particular are not quadratic families, although we must 
say that it is possible to get sometimes a generic point for some very specific and tiny 
sub-family of quadratic twists. 

3 This equality fixes once and for all our choice of the canonical height. Note that this 
height is twice the height in Silverman's book [Sil or in Krir's paper Kri so this explains 
the difference of a factor 2 between the formulae in this paper and theirs. 
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isomorphism between E and E d defined in the following way. Assume that 
the curves E and E d are given by: 

E : y 2 = X 3 + Ax 2 + Bx + c 
E d : y 2 = x 3 + Adx + Bd 2 x + Cd 3 



then ip d is: 



The non-trivial automorphism x i— ► x of Q(vd), which is the restriction of 
the complex conjugation if d < 0, acts by: 



4>d ■■ E £ d 

(x, y) i — ► (da:, d 3 / 2 y) 



MP) = -MP) (i) 

Whenever d and are coprime (and this will always be the case in our 
families), the conductor of E d is Nd 2 and we have: 

L(E d ,s) = ^2a(n)xd(n)n~ s 

n>l 

where Xd(-) — (") is the quadratic character associated to d. The sign of the 
functional equation satisfied by L(E d , s) is 

w(E d ) = w-Xd(-N). 

In the odd rank case (i.e. w(E d ) = —1), we are interested in the values at 
s = 1 of the derivatives of the L-functions. We have: 



n(E d ) c{E, 



E 



tors 



2 



where as usual Q(E d ) is the real period, R(E d ) is the regulator and c(E d ) = 
Yl p \Nd c p(P) ^ S the product of the local Tamagawa numbers. The Birch and 
Swinnerton-Dyer conjecture predicts that S(E d ) = \Ul(E d )\ if L'(E d , 1) ^ 
and S(E d ) = otherwise. 



2 Families of quadratic twists 

For each prime p dividing the conductor iV of E, we fix a sign tf p = ±1 so 
that ripiAr^p = w - We then define the set: 

T = jd < 0, fundamental discriminant with ^-^ = w p for all p \ iv| 
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and we let: 

F(T) = {def, \d\ < t} 

Then, our family of quadratic twists is the set (Ed) de:F and, for all these 
curves Ed, we have w(E^) = —1 by the above assumption on the product 
of the w p , s. It will be convenient for us to partition the family T into two 
subfamilies corresponding to the odd and even discriminant cases. Therefore 
we define: 

-^odd = yd G J 7 , d oddj and F dd{T) = |rf G FiT), d oddj 

Note that we will not need to consider the subfamilies corresponding to the 
even discriminants. 

For d G T with \d\ large enough, it follows from Proposition 2 of [De2j 
that, if we denote by c 4 the usual invariant of E (cf. |Cohlt §7.1]), we have: 



= s (2) 

d 8 (a,c 4 ) fijr lIcp(E d ) 

where 5 8 (d, c 4 ) = 2 if 8 | d and 2 | c 4 , and S 8 (d, c 4 ) = 1 otherwise, and 
fijc- is some positive number which does not depend on d. When L'(Ed, 1) 
is not zero then Ed(Q) has rank 1 and the regulator R(E<i) is equal to the 
canonical height h{Gd) of a generator Gd of Ed- So, the problem of studying 
the behavior of R(Ed) is roughly speaking the same as the one of studying 
the complexity of rational solutions of the associated Diophantine equations. 



2.1 On upper bounds for h(Gd) 

Lang's conjecture |Sill Conjecture 10.2] predicts that for a general elliptic 
curve E: 

i?(E)«|A min ( J E;)| 1 / 2+e 

where A min (E) is the minimal discriminant of E. In our family, we have 
A min (E d ) = <i 6 A min (i?) hence, this yields: 

R(E d ) « 

Of course, this upper bound is very far from what we really expect for our 
family. Indeed, using equation ^ and the fact that S(Ed) and c p (Ed) are 
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positive integers (so greater or equal to 1), the Lindelof hypothesis applied 
to L'(E d , 1) gives the following conditional upper bound: 

R(E d ) « e \d\^ 2+t 

In some cases, this upper bound can be proved on average. Anticipating 
on the results and notations of Section 3.1[ we prove: 



Proposition 1. Assume that N is square-free, L(E, 1) ^ and w p = +1 for 
all p | N. Then we have: 

' R(E d ) <T 1 / 2 logT (3) 



\F oM {T)\ ^ 

d ^ odd (T) 

Proof. This is a direct corol lary of a theorem of Ricotta and Vidick. Indeed, 



with the notations of section 3.1 



we have R(E d ) = h(G d ) < h(R d ) = Ah E (P L 



d) 



where He is the canonical height on E and P d G E[Qy d) is the Heegner point 



constructed in 3.1 Now, we apply the corollaire 3.2 of |Ri- Vi] . □ 



Remark. Classical conjectures predict that the number of discriminants d 
in our family for which L'(E d , 1) = should have density (we will come 
back to this fact later), so \J r odd (T)\ is roughly the number of terms in the 
sum of the formula above and hence the proposition really asserts that on 
average R{E d ) < \d\ 1/2+t for all d G T odd . 

2.2 On lower bound for R(Ed) 

Another conjecture of Lang asserts that h(G d ) ^> log |A min (.E d )|, thus we 
get: 

MG,)»log|d| (4) 
In fact, we have the more precise result: 

Proposition 2. If j(E) ^ 0, 1728, then there is an explicit constant M, 
depending on E and on the w p , such that we have for all d G T : 

h{G d )>^\og\d\ 
If w p = +1 for all p | iV, then one can take M = 1296 c(E) 2 . 
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Proof. We estimate \cm(c p (E d )) p \ Nd where c p (E d ) is the local Tamagawa 
number at the prime p dividing Nd. If p \ N, then c p (E d ) is either c p (E) 
if w p = +1, or c p (E*) if w p = — 1 where E* is any fixed twist of E by a 
discriminant that is not a square in Q p . If p | <i, then c p (E d ) is either 1, 2 or 
4. Hence, we have 

lcm^C^puv < 4 J] c p (£) J] c p (E*) 

p\N, w p =+l p\N,w p =-l 

Now, the result follows using Corollaire 2.2 of |Kri] and the fact that | A min (E d ) \ = 
|d| 6 |A min (£)|. □ 



Remark. 

1. With the same techniques, we can obtain similar results for j(E) = 
or 1728. 

2. One can prove (see for example |Sill exercise 8.17]) the following lower 
bound: 

h{G d )>±\og\d\+C (5) 

where C is some constant depending on E. The factor 1/3 in this 
formula is much better than the factor 1/M in Proposition [2j However, 
the constant C (which comes from the difference between the naive and 
the canonical heights) is negative and thus the estimate ^ is useless 
for small d (and in practice for all the gPs we are dealing with). On 
the other hand, the estimate of Proposition [2] is good enough for our 
applications and has no consequence on the main complexity of our 
method. 

3. The lower bound in Proposition [2] is optimal in the following sense: 
suppose that E is given by the equation y 2 = P(x) where P(x) is 
a degree 3 polynomial. Then, one can easily check that the point 
(rP(r), P{r) 2 ) belongs to Ep^(Q) and that the height of this point is 
w 4/31og|P(r)|. 

One expect much better lower bounds on average: it is proved in |De2j 
that predictions coming from Random Matrix Theory for derivatives of L- 
functions (see [Sna]) and Cohen-Lenstra type heuristics for Tate-Shafarevich 
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groups (see |Delj ) imply that for k > 0: 

J2 h ^ k » Tk,2 ~ e ( 6 ) 

L'(E d ,l)#0 

where the implied constant depends on E, k, e and w. 

2.3 Heuristics for the moments of R(Ed) 
For k > we let: 

L'(E d ,l)^0 

Equations Q and ([6]) imply that on average h(Gd) should be of the size of 
Ic?! 1 / 2 . In fact, one can make similar computations as in |De2j to estimate: 

£ R(E d ) k S(E d ) k 

der(T) 

L'(E d ,l)^0 

Then, Cohen-Lenstra type heuristics for Tate-Shafarevich groups (see |Delj ) 

predict that r^FTyJi S{Ed) k tends to a finite limit as T — ► oo whenever < 

k < 1. Therefore, using an empirical argument, we replace the term S(Ed) k 
by a constant and deduce the following heuristics: 

Heuristic for M k (T). For < k < 1 we have as T — > oo: 

M k (T) ~ A k T k ' 2 log^^+Wafc-i ^ 

for some constants A k and a k . 

The number a k comes from the contribution of the Tamagawa numbers in 
the Birch and Swinnerton-Dyer conjecture. More precisely we should have: 

• If E (or an isogenous curve) has full rational 2-torsion then a k = 4~ k . 

• HE has exactly one rational 2-torsion point (and no isogenous curve 
has full 2-torsion) then a k = \{4~ k + 2~ fe ). 
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For the other cases, we need to make the rather technical assumption that 
our restrictions on the discriminants are not incompatible with the use of the 
Chebotarev density theorem (see |De2j ). Then we should have: 

• If E has no rational 2-torsion point and its discriminant is not a square 
then a k = \A~ k + \2~ k + \. 

• If E has no rational 2-torsion point and its discriminant is a square 
then df, — | A~ k + |. 

Indeed, the equivalence §f§ depends only on the isogenous class of the curve, 
and this explains why we have to consider the curve in the class with the 
maximal rational 2-torsion point. 

If we restrict our family to negative prime discriminants, the effect of the 
Tamagawa numbers disappears and we have a k = 1. More precisely if we let: 

T' = |g? < 0,fund. disc, with f — J = w p for all p \ N and \d\ is prime | 

JF'(T) = {de r, \d\ < t} 

and 

d£T'(T) 
L'(E d ,l)jtO 

we expect the following heuristic: 

Heuristic for M' k (T). For < k < 1, we have as T — ► oo: 

M£(T) ~ A' k T k ' 2 log(T) fc ( fc+1 )/ 2 (8) 

Remark. These two heuristics are supported by our numerical data for the 
elliptic curves of conductor N < 17 as we will see in the last section. 

The asymptotics ([7| and ^ imply that on average the regulators of 
(E d ) de:F behave as ~ |<i| 1//2+e suggesting that 9 = e in the Saturday Night 
Conjecture (see [CRSWj ). From this we get a density of T 1 " 6 for the subset 
of d G F(T) such that L'(E d , 1) = 0, which is really surprising compared 
to the even-rank case. The numerical data seems to support this fact. On 
the other hand, extensive numerical computations by Watkins |Watj seem to 
indicate otherwise. Indeed we want to emphasize that one has always to be 
careful with deducing too strong of statements from numerical investigations. 
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3 Computation of generators 

We need to make a certain number of restrictions in order to be able to 
apply the method described in this section. First, we assume that E is the 
strong Weil curve in its isogeny class (in fact, we just need that the Manin's 
constant of E is equal to 1) and that j(E) ^ 0, 1728. These are just technical 
and not essential assumptions. Furthermore, we assume L(E, 1) ^ which 
implies that E(Q) has rank and that w — +1. This is a fundamental 
assumption and the method would not work without it. Finally, the family 
of discriminants T is obtained by taking w p = +1 for all p \ N. Hence, 
w(Ed) = — 1 and d is a square modulo 4iV for all d G T . 

The latter condition implies that one can apply the Heegner point con- 
struction to get a point Pd G E(Q(y/d)) of infinite order if L'(Ed, 1) ^ oQ For 
that one has to evaluate the modular parametrization at well chosen points 
r G X (N): 

ip : X (N) C/A -^E(C) 

n>l 

with X (N) = r (iV)\EI where r (A^) is the congruence subgroup of SX 2 (Z) 
of matrices with lower left entry divisible by N, HI = EI U Q is the completed 
upper half plane, A is the period lattice associated to E and p is the analytic 
isomorphism given by the Weierstrass function (and its derivative). 



3.1 Description of the method 

We now briefly describe the algorithm step by step. 

4 Classically the Heegner point method is used to construct directly a rational point on 
Ed(Q), see |Coh2[ Chapter 8.5]. However the direct construction of a point in a quadratic 
extension has been already done in connection with the problem of congruent numbers by 
N. Elkies, see [ElkJ. The main difference with the construction used in this article is that 
Elkics just wanted a strategy to compute efficiently a rational point of some quadratic 
twists of the elliptic curve 32a2, whereas we want to compute a generator of all the Ed(Q) 
for d £ T(T) of some large T. Hence, we really need to be careful in all the steps of the 
method in order to be as efficient as possible. We also have to use the full force of the 
Gross-Zagier formula and of the Birch and Swinnerton-Dyer conjecture in order to get as 
much information as possible all throughout our computations. 
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STEP 1. For each ideal class C in the class group Gl(d) ofQ(vd), we choose 
an integral ideal a € C such that: 

a = A Z + ~ B | — Z with N \ A and B = (3 (mod 2N) (9) 

where [3 = [3d is a fixed integer such that f3 2 = d (mod AN) . 
Then, to C = [a], we associate the Heegner point: 

-B + Vd 
T[a] = ~2A— 



Comments. The point T[ a ] lies in the upper half plane and is a well defined 
point in X (N). Nevertheless, in order to make the computations as easy 
as possible, we need to choose a such that A is as small as possible. Us- 
ing classical algorithms (see |Cohlj ). we can compute a set of ideals {di}i 
representing all the classes of C\(d): 



—hi + yd 
dj = ctjZ H Z 

with < ai <^ \d\ l l 2 where the implied constant is explicit. We can assume 
without loss of generality that the a^s are relatively prime with N. Then, 
the ideals ajtt satisfy Q where 

n = iVZ+^^Z 
2 

From this it follows that one can choose the ideals Oi's in such a way that we 
have the following lower bound: 

3(r w ) » l/N (10) 

The complexity of this step is thus dominated by the class number of Q(y/d), 
hence is at most 0(|ei| 1//2 log |c?|). 

STEP 2. We compute 

where the sum is over the classes of Cl(d), and then a complex approxi- 
mation of Pd = p{zd) G E(C). The theory of complex multiplication and 
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of Heegner points imply that Pd 6 E(Q(yd)). Using this approximation, 
we try to recognize the four rational numbers r\, s\, r 2 and S2 such that 
Pd = {ti + siVd, r 2 + s 2 \fd) and test if Pd is a point of infinite order. 

Comments. This is the main step of the method. Note that one can reduce 
the number of evaluations of by 2 using the following trick. Once we have 
already computed ip(r^), since w = +1 we can deduce from it <£>(T| -i n ]) 
using the formula: 

f(r [a] ) = -<p(T [a -i n] ) + Q (11) 

where Q is an explicit rational torsion point in E(Q) depending only on E. 

Given a complex number xp d that is an approximation of the x-coordinate 
xp d of the point Pd computed as explained above, we need to recover from 
it the two rational numbers r± and si such that xp d = r\ + swfd. Note that 
for candidate values r x and Si, one can check if they are indeed correct by 
trying to compute two rationals r 2 and s 2 such that (r x + Si \[d, r 2 + s 2 \fd) G 
E(Q(Vd)). Let f = 9ft(x P J and s = ^s(x p d )/ \fd. For e > 1 we look for a small 
integral relation (using the LLL-algorithm) between the columns Ci, 62,63 
of the matrix 

-10 e L 1Qe ^l 
-io e LiO's] 

1 

where [•] denotes the closest integer. Indeed, for such a relation, say 

\\G\ + A2C2 + A3C3 

of norm M, we have that A1/A3, resp. A 2 /A 3 , is an approximation of f, 
resp. s, with an error less than v^f/10 6 , and the denominator A3 is smaller 
(in absolute value) than \fM . In order for this method to work, we need 
to compute f and s at a suitably large enough precision and to choose e 
accordingly. More precisely, to recognize xp d as an element of Q(\/d) we 
need about h(P d ) di gits. Bounding the coefficients a{n)/n by 1 in the sum 



defining and using (10), we see that we need to sum approximatively h(P d ) 
coefficients for (p. The Gross-Zagier theorem Gro-Zag] asserts thatj^] 



h{Pd) ~ I^oKE) (12) 

5 Actually, the Gross-Zagier theorem only applies for odd d's. For even d's the formula 
is a conjecture of Hayashi [Hay]. 
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Applying the Lindelof hypothesis we deduce that h(Pd) <C |<i| 1//2+e . Hence, 
the complexity of this step is <C | cZ| 1 / 2 -+- e | Cl(c?) | <C This step can fail in 

two ways. First case: the computation has not been done to a large enough 
precision. In that case we have to increase the precision and start over. 
Second case: the point P d is a torsion point and in that case L'(E d , 1) = 0. 
If we suspect Pd to be in fact a torsion point, we can compute directly 
an approximation of L'(Ed, 1) and prove that it is indeed zero using the 
following proposition (whose proof we postpone to after the proof of the 
next proposition). 

Proposition 3. If 

then L'(E d , 1) = 0. 



STEP 3. If Pd is a point of infinite order, i.e. STEP 2 has succeeded, then 
the point Rd = ip(Pd — Pd) is a point of infinite order in E d (Q). We divide 
it in the Mordell-Weil group E(Q) until we get a generator Gd € E d (Q). We 
define the integer t d by R d = t d G d (mod E d (Q) tor ). 



Comments. The point Rd is rational since R d = 4>(Pd) + ip{Pd) by (JB. If 
L'(E d , 1) 7^ then we know that G d is a generator of Ed(Q) modulo torsion. 

Proposition 4. h(Rd) = AfiE{Pd), hence Rd is non-torsion if and only if Pd 
is non-torsion (that is if and only if L'(Ed, 1) ^ 0). 

Proof. The height does not depend on the model of the elliptic curve, hence 
h(Rd) = hE^Pd — Pd)- Furthermore, equation (11) implies that P d = —Pd 



plus a rational torsion point. □ 
Proof of proposition^ We use the lower bound from proposition [2] for h(R d ) 



and equation (12). □ 



From Proposition [2] we know that: 



h(R a 



K,/|<36r(A-)J^|<|d| 1 /^ (13) 

Hence there are finitely many primes p for which we need to check p-di- 
visibility. Also, it is well-known that Ed(Q)tors does not depend upon d (for 
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all cfs except at most one) and can only be ~ {0}, Z/2Z or Z/2Z x Z/2Z. 
Therefore we need to be careful about torsion only when we consider 2- 
divisibility which can be tested easily using 2-division polynomial. For an 
odd prime p, we use the following method to rule out p-divisibility. We find 
a prime r, of good reduction, such that the order a of the group E d (F r ) is 
divisible by p. Then if (a/p)R d is not zero in E d (¥ r ), we know that R d is not 
divisible by p in E(¥ r ), and thus in E d (Q) too. If after having performed a 
large number of such tests, we have not been able to prove that Rd is not 
divisible by p, then we "know" that the point must be divisible by p and we 
perform the division]^] 

STEP 4. We compute the regulator of E d (in the rank 1 case) which is equal 
to h(G d ) and the order of the Tate-Shafarevich group \m(E d )\. 

Comments. We can compute the order of using: 

Proposition 5. Under the Birch and Swinnerton-Dyer conjectur^\ the fol- 
lowing equality holds 

\m(E d 

lu \ il 1 1 cyn, ) - i s(j \ iA min yn ))o 8 ya,c 4 ) 

P \d 

where sg(x) = 1 if x < and sg(x) = 2 otherwise. 
Proof. Indeed, we have: 



\E(Q) tor \ 


2 


E d (Q) tor \ 


2 p2 




m(E)\c(E) 2 


2sg(A mm (E))5s(d,c 4 )Uc P (E d ) 



Fh(C ) - MP 1 - HEA)L'(E d ,l)y/\d\ 
t d h{G d ) - Ah(P d ) — ^ 



Now we replace L(E, 1) and L'(E d , 1) by the values predicted by the Birch 
and Swinnerton-Dyer conjecture. After simplifying the regulator h(G d ) on 
both sides, we get: 



ME d )\ 


\m(E)\c(E) 


\E(Q) tov \ 


2 


E d (Q)tor\ 


2 



voiyn/ ) 



Since for all p \ N we have w p = +1, the curves E and E d are isomorphic 
over Q p and thus c p (E d ) = c p (E). So c(E d ) = c(E) rip|d c p(^)- Finally a 



6 Indeed, in all cases, either we could prove that the point is not divisible by p by such 
a test, or we could actually divide it by p. 

7 For even c?'s, we need again to assume the conjecture of Hayashi [Hay]. 
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computation of the periods of E d shows that: 

n(E) n(E d ) y^T = r 2 s 8 (d, C4 ) if A min (E) < 

vol(E) \ Ss(d, Ci ) if A min (£) > 

and the proposition follows. □ 

Remark. The order of the Tate-Shafarevich group is a square, therefore the 
proposition implies that the following quantity must be a square: 

2 sg(A min (E)) 5 8 (d, c 4 ) JJ c p (E d ) 

P \d 

From the above we see that for each individual d the complexity for 
computing h(G d ) and | HI ( 2^)1 is at worst 0(\d\ 1+e ). From these values we 
can deduce the value of L'(E d , 1) at arbitrary precision. Note that the direct 
computation of L'(E d , 1) by the rapidly converging series needs also 0(|d|) 
termsj^] Nevertheless, for large precisions, in practice, it is often much more 
efficient to compute L'(E d , 1) as a by product of our computations than to 
evaluate it directly. This is probably due to the fact (see the discussion 
on the computations) that the implied constant is small in the prediction 
Ml(T) = 0(T 1 / 2 (logT) a ). 



3.2 An example 

We take E = Hal : y 2 + y = x 3 - x 2 - Wx - 20 and d = -79 so that the 
curve E d has minimal equation: 

E d : y 2 + y = x 3 + x 2 - 64490x + 11396008 

We take (3 = 3 so that (3 2 = -79 (mod 44). The class group Cl(-79) of 
—79) is cyclic of order 5, and the ideals: 

-3 + v^79 -3 + v^79 
a=HZH Z, b = 22Z H Z, 

_3 _|_ ^/_70 

c = 44Z+ Z, a- x n and fr x n 



8 More generally, in order to compute L'(E, 1) for an elliptic curve E, one needs to sum 
the first 0(y/N) terms of the series, where N is the conductor of E, and the constant in 
the "O" depends on the required accuracy. 
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where 

n = iVZ H Z = 11Z H Z 

2 2 

form a complete set of representatives of the ideal class group. We compute 
z = 23?(0(r [a] ) + 0(r [b] )) + 0(r [c] ) G C/A 

and we find 

p d = p( z ) « (-3.59000 • • -+0.22000 • ■ • V /Z 79, 5.17600 • • -+0.61600 ■ ■ ■ V^79) 
so we easily recognize 

Pt = (iH!^, 647 - *(Q<^)) 



From this, we get the point R d = ip(P d ) + ^(P d ) = (47,2910) G £ d 
And Formula (13) says that \£ d \ < 293 where Gd = C-dRd- We find that 
the point Rd is divisible by 2, more precisely Rd = —2(26,3120), so that 
(26,3120) = l' d G d with \Q < 73. We then easily check that the point 
(26, 3120) is not divisible by any prime < 73 in the group Ed(Q), hence one 
can take Gd = (26,3120) and \£ d \ = 2. Proposition ^ gives: 

ME d )\ = 1 



4 Discussion and numerical data 

We have computed, using the method described in the previous section, the 
regulators and the order of the Tate-Shafarevich groups of the twists E d of 
E of the four elliptic curves llal, 14al, 15al and 17al, and for all available 
discriminants d G JF(1.5 x 10 6 ) with w p = +1 for all p | N. We discuss in 
this section the data we obtained and compare it with the heuristics. All the 
computations have been performed using the PARI/GP system [PARI] and 
the data is available at 

http : //math.univ-lyonl . f r/~roblot/tables .html 

We begin with the curves of prime conductor (llal and 17al) since for 
the last two curves (14al and 15al), the congruence conditions are more 
restrictive and therefore the number of discriminants in JF(1.5 x 10 6 ) is quite 
small compared to 1.5 x 10 6 . 
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4.1 The curve llal 

The curve E is defined by E : y 2 + y = x 3 — x 2 — lOx — 20. It has conductor 
N — 11 and rank over Q. We have w\\ = +1. 



4.1.1 Numerical results for all discriminants 

• Number of discriminants: |J^(1.5 x 10 6 )| = 208977. 

• Largest regulator: « 9945 (for d = -1482139). 

• Number of extra- vanishing: 638. 

We have E(Q) tors ~ Z/5Z and there is no curve in its isogeny class having 
rational 2-torsion. Hence the heuristics predict that: 

M fc (T) ~ A fc T fc/2 log(T) M ^ + MS + 5^-i 

for some constant A\~. We computed A^ numerically to fit the data (values 
found: Ay a, « 0.60, A 1/2 w 0.33, A 3/4 « 0.16, Ai w 0.07) and we plot the 
graph of the function given by the heuristics and the points (T, Mfc(T)) for 
T = 1, 2, . . . , 150 x 3 • 10 4 and for k = 1/4, 1/2, 3/4 and 1. As it can been 
seen the graphs (see Figure [2]) are in close agreement. 



4.1.2 Numerical results for prime discriminants 

• Number of prime discriminants: 28535. 

• Largest regulator: ps 9250 (for d = —1433539). 

• Number of extra-vanishing: 0j^] 

The heuristics for prime discriminants predict that: 

M' k (T) ~ A , k T k ' 2 \og(T) k{k+1 V 2 

for some constant A' k . We computed A' k numerically to fit the data (values 
found: A'^ 2 ~ 0.20, A[ ps 0.03) and we plot the graph of the function given 
by the heuristics and the points (T, M' k (T)) for T = 1, 2, ... 150 x 3 ■ 10 4 , and 
k = 1/2, 1 (see Figure g. 

9 Thcre is no extra- vanishing in this case using the results of |An-Bu-Fr] . 
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GflQOQO 80O00O le+06 1.2e+06 1.4e+06 

I Discr inlnants I 




le+06 1.2e+66 1. le+06 



I Discriminants I 



Figure 1: Regulators of the twists of llal. On the left, all the regulators 
are plotted, and on the right, only regulators less than 10. The gap between 
the x-axis and the minimal heights is clearly visible on the right. 
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Figure 2: Moments of order i, |, | and 1 of the regulators of the twists of 
llal and the function given by the heuristics. 
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vanishing. data" 
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ordre 1/2 X 

ordre 3/4 * 
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□ 

□ 
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□ 

1.2 - ****** 
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I Discrihinants I 



L'(E d , 1) for E = Hal. 



Figure 3: Number of extra- vanishing of Figure 4: Moments of order |, |, | and 1 

for the order of the Tate-Shafarevich groups 
of the twists of Hal. The heuristics suggest 
that the moments of order k < 1 tend to a 
constant (depending on k) whereas the mo- 
ment of order 1 should tend to infinity. 




Figure 5: Moments of order | and 1 for the regulators of the twists of Hal 
by prime discriminants and the functions given by the heuristics. 
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4.2 The curve 17al 

The curve E is defined by E : y 2 + xy + y = x 3 — x 2 — x — 14. It has 
conductor 17 and rank over Q. We have wn = +1. 

4.2.1 Numerical results for all discriminants 

• Number of discriminants: 215305. 

• Largest regulator: « 31746 (for d = —1257787). 

• Number of extra-vanishing: 1140. 

Remark. Note that the graphs of extra-vanishing for the curves llal (Fig- 
ure^ and 17al (Figure [8]) suggest that the density of extra- vanishing is larger 
for the twists of 17al than for those of llal. However the asymptotic for 
the moments of the regulators is smaller (as T — > oo) for 17al than for llal 
which suggest that there are more constraints on the regulators of the twists 
of llal and thus imply in turn that we should have more extra- vanishing 
for this family. In fact, the constants Ak in the asymptotics of Mfc(T) are 
larger for the curve 17al, but asymptotics of the functions Mf.(T) for the 
curve llal are larger than for the curve 17al for very large values of T that 
are completely out of reach for computations. Therefore our guess is that 
the density of extra- vanishing for the twists of llal will become greater than 
that for the twists of 17al for those very large values. 

The curve 17a2 has full rational 2-torsion, hence the heuristics predict 
that 

i /ri k(fc+i) i . 

M k (T) ~ A k T k ' 2 \og(T)— + ^- 1 

for some constant A k . We computed A k numerically to fit the data (values 
found: A 1/A w 0.97, A 1/2 « 0.75, A 3/4 w 0.47, A 1 w 0.25 and we plot the 
graph of the function given by the heuristics and the points (T,M k (T)) for 
T = 1, 2, ... 150 x 3 • 10 4 , and k = 1/4, 1/2, 3/4 and 1 (see Figure Q. 

4.2.2 Numerical results for prime discriminants 

• Number of prime discriminants: 28601. 

• Largest regulator: ps 31745 (for d = —1257787). 
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Figure 6: Regulators of the twists of Hal. On the left, all regulators are 
plotted, and on the right only regulators less than 10. The gap between the 
and the minimal heights is cleary visible on the right. 




Figure 7: Moments of order j, § and 1 of the regulators of the twists of 
17al and the functions given by the heuristics. 
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Figure 8: Number of extra- vanishing of 
L'(E d , 1) for E = Hal. 




Figure 9: Moments of order |, |, | and 1 
for the order of the Tate-Shafarevich groups 
of the twists of 17al. The heuristics suggest 
that the moments of order k < 1 tend to a 
constant (depending on k) whereas the mo- 
ment of order 1 should tend to infinity. 




Figure 10: Moments of order | and 1 for the regulators of the twists of 17al 
by prime discriminants and the functions given by the heuristics. 
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• Number of extra-vanishing: 0, 

The heuristics for prime discriminants predicts that: 

M' k (T) ~ A' k T k/2 \og(T) k{k+1)/2 

for some constant A' k . We computed A' k numerically to fit the data (values 
found: A'y 2 ~ 0.41, A[ « 0.12) and we plot the graph of the function given 
by the heuristic and the points (T, M' k (T)) for T = 1, 2, ... 150 x 3 • 10 4 , and 



k = 1/2, 1 (see Figure 10). 



4.3 The curve 14al 

The curve is defined by E : y 2 + xy + y = x 3 + Ax — 6. It has conductor 
iV = 14 and rank over Q. We have w 2 = w-j = +1. 



• Number of discriminants: 66516. 

• Largest regulator: « 16937 (for d = —1416631). 

• Number of extra-vanishing: 262. 



We have E(Q) toIS ~ Z/3Z, and there is no curve in the isogeny class 
having full rational 2-torsion. Hence the heuristics predict that 

M k (T) ~ ^T^log^^+I^^ 2 ")- 1 



for some constant A k . We computed A k numerically to fit the data (values 
found: Ay a « 0.82, Ay 2 w 0.56, A 3/4 w 0.33, A x w 0.17) and we plot the 
graph of the function given by the heuristics and the points (T, Mfc(T)) for 
T = 1, 2, ... 150 x 3 ■ 10 4 , and fc = 1/4, 1/2, 3/4 and 1 (see Figure pk. 



4.4 The curve 15al 

The curve E is defined by E : y 2 + xy + y = x 3 + x 2 — lOx — 10. It has 
conductor iV = 15 and rank over Q. We have W3 = W5 = +1. 

• Number of discriminants: 71254. 

10 Thcrc is no extra- vanishing in this case using the results of |An-Bu-Fr] . 
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Figure 11: Number of extra- vanishing of 
L'(E d , 1) for E = Ual. 



Figure 12: Moments of order |, |, | and 1 
of the order of the Tate-Shafarevich groups 
of the twists of 14al. The heuristics suggest 
that the moments of order k < 1 tend to a 
constant (depending on k) whereas the mo- 
ment of order 1 should tend to infinity. 
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290000 480000 600000 800000 le+OG l,2e*BG 1.4e+06 l.Ge+O 

Figure 14: Number of extra-vanishing for 
L'(E d , 1) for E = 15al. 



Figure 15: Moments of order z, h> f an d 1 
of the order of the Tate-Shafarevich groups 
of the twists of 15al. The heuristics suggest 
that the moments of order k < 1 tend to a 
constant (depending on k) whereas the mo- 
ment of order 1 should tend to oo. 





200000 400000 GO0OQO 000000 le+OG 1.2e+86 1.4B+06 1.6e+0 200000 400000 GO0OQO 000000 le+OG 1.2e+86 1.4B+06 l.Ge+0 

Figure 16: Moments of order |, |, | and 1 for the regulators of the twists 
of 15al and the functions given by the heuristics. 
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• Largest generator: pa 19352 (for d = —1297619). 

• Number of extra-vanishing: 406. 

We have E(Q) toIS ~ Z/4Z x Z/2Z, hence it has full 2-torsion. The 
heuristics predict that 



M k (T)~A k T k / 2 \og(T) 



fc(fc+l)/2+4- 



for some constant A k . We computed A k numerically to fit the data (values 
found: A 1/4 pa 0.97, A 1/2 pa 0.75, A 3/4 pa 0.47, A 1 pa 0.25) and we plot the 
graph of the function given by the heuristic and the points (T, Mfc(T)) for 



1, 2, ... 150 x 3 ■ 10 4 , and k = 1/4, 1/2, 3/4 and 1 (see Figure 16) 
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